Crossed Andreev reflection in a d-wave superconductor 
with two quantum point contacts 
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We theoretically study the crossed Andreev reflection in a hybrid nanostructure which comprises a 
d-wave superconductor and two normal-metal quantum wires. When the superconductor of the (110) 
oriented surface is in contact with the wires parallel and placed close to each other, the Andreev 
bound state is formed by the crossed Andreev reflection. When the contact barrier potential is 
sufficiently large, two sharp peaks appear in the conductance well below the gap structure, which 
originate from the bonding and antibonding Andreev bound states. We propose that these Andreev 
bound states form a two- level quantum system (qubit). 

PACS numbers: 74.45.+c,74.50.+r,81.07.Lk,73.40.Gk 
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Quantum transport in nanostructures is of current in- 
terest in both theoretical and experimental studies. In a 
quantum wire, its width is so narrow that electron waves 
are strongly confined in the transverse direction and their 
transverse momentums are quantized. In a ballistic con- 
duction of electrons through a quantum point contact 
with width comparable to the Fermi wavelength, a step- 
like structure with step of 2e 2 jh appears in the conduc- 
tance as a function of Fermi energy or width [1-3] . 

One of the fundamental consequences of superconduc- 
tivity is the Andreev reflection at the interface of a 
normal-metal and a superconductor (SC) [4,5]. This phe- 
nomenon corresponds to an incoming electron from the 
normal side being reflected as a hole, thereby adding a 
Cooper pair in the superconducting condensate. In a 
tunnel junction of a normal-metal and a (110) oriented 
rf-wave SC, the zero bias conductance peak appears due 
to the formation of the Andreev bound state at the in- 
terface [6,7]. However, when a single quantum wire of a 
single conducting channel is in contact with d-wave SC 
of the (110) oriented surface, the Andreev reflection is 
completely suppressed due to the quantum mechanical 
diffraction of electron waves at the narrow opening [8]. 

A basic question arises what happens if two quantum 
wires are in contact to the (110) oriented d-wave SC (see 
Fig. 1). When an electron is injected into SC from one 
of the wires, there is a possibility that the Andreev hole 
is reflected back into another wire due to the non-local 
effect called the crossed Andreev reflection (CAR) [9,10], 
thus providing an ideal system to study CAR. 

In this Letter, we explore quantum-interference effects 
due to the crossed Andreev reflection in a d-wave SC with 
two quantum wires that are parallel and placed close to 
each other. It is shown that the resonance peak in the 
conductance is split into two sharp peaks at low energies 
when the barrier potential of the contact is sufficiently 
large. The lower and higher energy peaks correspond to 
the bonding and the antibonding Andreev bound states, 
respectively. This suggests that these Andreev levels 



form a two-level quantum system (qubit), whose pop- 
ulation can be controlled by application of bias voltage 
and/or electromagnetic field, exhibiting a coherent oscil- 
lation of the Andreev qubit (Rabi oscillation). 

We examine the quantum transport in a hybrid nanos- 
tructure of a d-wave SC and two normal-conducting 
quantum wires. Figure 1 shows a model structure in the 
x-y plane with a two-dimensional (2D) d-wave SC occu- 
pying the left half space, and two quantum wires of width 
w, lead 1 and lead 2, which are parallel along x and con- 
nected to SC at y = ±L/2. The wave functions of elec- 
tron and hole like quasiparticles with excitation energy 
E in the electrodes are determined by the Bogoliubov-de 
Gennes equation 
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where Hq = — ?TV 2 /2m — ep is the single-particle Hamil- 
tonian with the Fermi energy €f = h 2 kp/2m. For sim- 
plicity, the Fermi wave number fcp and the effective mass 
m are common for all electrodes, and the amplitude of 
the gap function A(k, r) is uniform in SC and vanishes 
in the wires. In an anisotropic SC, the pair potential 
A(k, r) is a function of wave vector k, and its value is 
defined on the Fermi surface in the direction of k. In 
leads 1 and 2 with infinite wall boundaries, the wave 









Lead 1 « •— *■ 




















L 

( + 1 




Lead 2 -*-^>h 










a!-wave superconductor 



FIG. 1. Schematic diagram of two- leg normal-conducting 
quantum wires in contact with a d-wave superconductor. 
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function of an electron propagating along the x axis with 
wave number k in the nth electron subband is given by 
the product of the plane wave e lkx and the nth transverse 
wave, Xn(y — L/2) in lead 1 and Xn(y + L/2) in lead 2, 
where Xn(y) = (2/w) 1 / 2 sin [(nir/w) (y + w/2)], and has 
the energy E = (H 2 /2m)[k 2 + (mr/w) 2 } - e F . In the 
following, we restrict our calculation to the case where 
only the lowest subband (n = 1) is occupied by electrons 
(or holes) in the wires, which is realized for those wires 
satisfying 7r < k-pw < 2tt. When an electron with energy 
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incident from lead 1 into SC, the wave functions in leads 
1 and 2 are given by ^\(x,y) = (pi(x)xi(y — L/2) and 
$2(x,y) = f 2 {x)xi(y + L/2), where 
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Here, r\\ and r\ \ are the amplitudes of the normal re- 
flection (NR) and Andreev reflection (AR), respectively, 
while and rf£ are those of the crossed normal reflec- 
tion (CNR) and the crossed Andreev reflection (CAR), 
respectively. A similar treatment is made for an inci- 
dent electron from lead 2. Since E < A <C e F , we put 
ki w fcp[l — (Tr/fcpw) 2 ] 1 / 2 in the following. 

We employ the Blonder-Tinkham-Klapwijk (BTK) ap- 
proach [5] to calculate the conductance of the structure. 
To carry out the calculation analytically, we make the 
Andreev approximation that neglects all the evanescent 
modes [11], and put the wave function of SC in the form 

Hf a (x,y) = |" F tf( Py ) ( r X + ) eV^V^dp,, 

+ f[ tf(Py) (^f) e-'VtF^'jPWdpy, (4) 

where T+ = A* + /(E + fl + ) and T_ = A_/(£ + Q_) 
with fi± = ^JE 2 - |A±| 2 [6-8]. The first and sec- 
ond terms in Eq. (4) are the transmitted QP waves 
on the electron-like and hole- like branches, respectively. 
In SC of d-wave symmetry, QPs in different branches 
feel different pair potentials, A + = Ao cos 2(6 — a) and 
A_ = Ao cos 2(9 + a), where a is the angle between the 
(100) axis of SC and the normal to the interface (see 
Fig. 1), and 6 = sin -1 (p y /k F ) is the propagation angle 
relative to the x axis. 

In the following, we focus on the (110) oriented surface 
of SC, i.e., a = 7r/4, as shown in Fig. 1. In this case, QPs 
in the electron-like and hole-like branches move in the 
pair potentials of opposite sign, A± = ±A, where A = 
2A (p v /k F )^l - (p y /k F ) 2 , and ft± = ft = VE 2 - A 2 . 
The barrier potential at the interface between the wires 
and SC is taken into account by the 5-function-type po- 
tential with amplitude (h 2 k F /2m)Z , Z being a dimen- 
sionless parameter [5]. The boundary conditions for the 



wave functions at the interfaces are *£ s (0,y) = ^i(0,y) 
and [d x ^ s (x,y)-d x ^i(x,y)} x=0 = k F Z^i(0,y) (i = 1,2) 
appropriate for the 5(x) potential. The matching tech- 
nique to the boundary conditions [3] yields the reflection 
coefficients 
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with fci = fci/fcp, V± = (k 1 +T) 2 - (g c ±iZ) 2 -g 2 , and 

r = f[ |§Vi-b/fc F )Vb), 0) 



dp Q, 
2^E 



^l-(p/k F ) 2 <p 2 (p)cos(pL), (10) 



g s = i 



l-(p/k F ) 2 <p 2 (p)sm(pL), (11) 



where ip(p) = (p\xi) is the overlap integral of xi{y) an d 
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Note that T in Eq. (9) represents the local coupling and 
is independent of distance L, while g c and g s represent 
the non-local coupling and are dependent on L. In the 
limit of L — > oo, where the two contacts are independent 
(g c = g s = 0), one has rff = (k 1 -T-iZ)/(k 1 +T+iZ), 
r eh _ r ee _ r eh _ recover i n g the complete suppression 
of AR in a single quantum wire with a single transverse 
mode [8]. The reflection coefficients r%2, r^, r|f , and rfi 
for an incident electron from lead 2 are obtained from 
those in Eqs. (5)-(8) by the replacement 1 <-> 2 and L ~ > 
— L. It follows from Eqs. (6)-(8) that, when there is no 
barrier potential (Z = 0) at the interface, AR is absent 
(rf{ = 0), whereas CNR and CAR are finite with the 
ratio r\l/r<&=g c /g s . 

When bias voltage V is applied to the two leads, the 
conductance G at zero temperature (T = 0) is given by 



4e 2 



Vll\ 2 



r ee \ 2 
r 12l 



„eh\2 



\rl h 2 \ 2 ) E = 



(13) 



where h is the Planck constant. Figure 2 shows the con- 
ductance G vs V for k F w = 4, k F L = 8, and different 
values of Z. The conductance is normalized by the nor- 
mal state value Gn — G(A = 0), which takes Gn ~ 
16(e 2 /h)(k 1 /k F )T N /Z 2 (Z > 1) with T N = T(A = 0). 
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FIG. 2. Normalized conductance G/Gn as a function of 
bias voltage V for different values of interfacial barrier pa- 
rameter Z. Inset shows G/Gn vs normalized voltage ZV . 
The circles in the inset are calculated from Eq. (16). 

For a low barrier potential of small Z, the conductance 
decreases monotonically with decreasing eV below Ao. 
As the barrier potential becomes higher, a peak struc- 
ture appears well below Ao and shifts towards lower eV, 
developing the double peak structure with increasing Z. 
If the conductance is plotted as a function of normalized 
voltage ZV as shown in the inset, the conductance peaks 
fall into the same position, indicating that the resonance 
peak positions are scaled by 1/Z. 

Let us examine the origin of the double peak struc- 
ture in the conductance for the tunneling case (Z ^> 1). 
Since w i\A\ for E <C A , T , Q c , and Q s in Eqs. (9)- 
(11) have the forms: T w ifA /E, Q c w ig c A /E, 
and Q s = ig s A /E, where / = {ET '/iA ) E ^o, g c = 
(Eg c /iA ) E ^o, and g s = (EQ s /iA ) are energy inde- 
pendent quantities, so that the Andreev reflection coeffi- 
cients are calculated in the resonance forms 
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where g = g,/2f, E ± = (/ ± \g c \ + «? 2 /2/)A /Z, and 
7 = (ki/kp)fAo/Z 2 . The other coefficients are r\\ w 
— 1 — gr^ and rf| ~ 0(g 2 ). The resonance energies -E± 
and intensity g depend strongly on lead separation L, ex- 
hibiting damped oscillations with the period of the Fermi 
wave length Af = 27r/fcp. We note that the resonance po- 
sitions E± and their separation E + — E- are scaled by 
1/Z, while the line width 7 is scaled by 1/Z 2 . For Z > 1, 
the line width of the peaks is much smaller than the sepa- 
ration, and therefore a well-separated two-peak structure 
is formed as shown in the inset of Fig. 2. In this case, the 
conductance G is written as the sum of two Lorentzians 
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where X = (irg 2 f /4^v). The simple formula (16) repro- 
duces the numerical result in the inset of Fig. 2, if the 
calculated values (/ = 0.405, g c = -0.036, g s = -0.0345, 
and Tn = 0.68) are used in Eq. (16). Note that the peak 
height of G/Gn increases in proportion to Z 2 . 

To elucidate the formation of the Andreev bound 
states, we calculate the QP wave function ^ s (x,y) — 
in SC, where \£f and are the electron and 
hole wave functions, respectively. Figure 3 shows the 
mapping of the absolute squares, l^l 2 and \^' s l \ 2 , at the 
resonance energies on the xy plane for kpw = 5, k-pL = 8, 
and Z = 50, when an electron is incident from lead 1. It 
is clearly seen that the QP wave functions are strongly 
localized with very large peaks (red color) near the con- 
tacts due to the formation of the Andreev bound states, 
and that the decaying QP waves are traced to the (010) 
and (010) directions along which multiple reflections of 
electron and hole take place. It is also seen that the 
formation of the localized states, especially the relative 
position of the peaks with respect to the contacts, is dif- 
ferent between the electron and hole QPs, and between 
the resonance energies E±. For Z ^> 1 and around E±, 
the explicit form of ^ s inside SC is obtained as 
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where, except very close to the interface (x rs 0), 
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with q — A/fcp — p 2 and A = 2Aopq/k F . At lower res- 
onance energy £L, \t s is dominated by the first term 
whose electron (hole) wave function -ip- (ip—) is an odd 
(even) function of y. At higher resonance energy E + , ^ s 
is dominated by the second term whose electron (hole) 
wave function tp^_ (iplL) is an even (odd) function of y. 
These results indicate that the electron (hole) wave func- 
tions at the lower and higher bound states have different 
parity with respect to y. It is noteworthy that the An- 
dreev holes are reflected back into leads 1 and 2 in phase 
at £L and out of phase at E + (cf. Eqs. (14) and (15)). 

The Andreev bound states appeared in the present sys- 
tem have the following implication. When the interface 
barrier potential is very high (Z 3> 1), the current flows 
through the system via the Andreev bound states with 
very sharp energy levels at E±. When the bias voltage 
V is set between V- = E- /e and V+ = E + /e of the con- 
ductance peaks, the lower bound state at E_ is occupied 
by Andreev quasiparticles while the higher-energy bound 
state at E + is empty; the lower bound state is viewed as 
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FIG. 3. Mapping of Andreev bound states on the xy plane, 
when an electron is incident from lead 1 to a d-wave SC. (a) 
and (b) are the absolute square of electron wave functions 
| VEf J | 2 at energies E- and E+; (c) and (d) are that of hole 
wave functions | 2 at E- and E + . 

the ground state and the upper bound state as the ex- 
cited state, and therefore the system in this setup forms 
a two- level quantum system (qubit). For simplicity, we 
denote the ground state and excited state by |0) and 
|1), respectively. The Andreev bound states of electron 
(hole) QP at E- and E + have dipole moments (l|y|0) 
due to the different parity of the states with respect to 
y. Therefore application of the electric field eE cos(wi) 
oscillating along y causes a significant Rabi oscillation 
with frequency V = eEo(l\y\0) /h between the Andreev 
bound states, if u is tuned to the resonance condition 
u> = E + — E- and the damping rate 7 of the states is sub- 
stantially smaller than TiV . The significant coherence os- 
cillation requires a smaller damping 7, which is achieved 
by inserting tunnel barriers at the interfaces. The con- 
trolled evolution between the two states |0) and |1) is 
realized by applying resonant microwaves to the system. 
The two level Andreev bound states in the proposed sys- 
tem provides a new possibility of superconducting qubit 
(d-wave Andreev level qubit) [12]. 

The realization of quantum two- level systems (qubit) is 
fundamental issues in both physics and information tech- 
nologies, because qubits are basic elements for quantum 
computation [13]. Recently, the coherent manipulation 
of the states in two-level systems has been demonstrated 
in nanostructured superconducting circuits [14,15]. 

In summary, we have theoretically studied the 
quantum-interference effects caused by the crossed An- 
dreev reflection in a hybrid nanostructure which com- 



prises a (i-wave superconductor and normal-conducting 
quantum wires. When a superconductor of the (110) 
oriented surface is in contact with the two-leg quantum 
wires via tunnel barriers, the resonance bound states are 
formed at low energies due to the crossed Andreev reflec- 
tion. As a consequence, two well-separated sharp peaks 
appear in the conductance well below the superconduct- 
ing gap structure. The lower and higher conductance 
peaks correspond to the bonding and the antibonding 
Andreev bound states whose wave functions have differ- 
ent parity. We propose that these Andreev bound states 
form a two-level quantum system (d-wave Andreev level 
qubit), which can be controlled by application of bias 
voltage and/or radiation of oscillating electric field. 
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